Study of transients in the propagation of nonlinear waves 
in some reaction diffusion systems 
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We study the transient dynamics of single species reaction diffusion systems whose reaction terms 
f(u) vary nonlinearly near u « 0, specifically as f(u) « u 2 and f(u) w u 3 . We consider three cases, 
calculate their traveling wave fronts and speeds analytically and solve the equations numerically 
with different initial conditions to study the approach to the asymptotic front shape and speed. 
Observed time evolution is found to be quite sensitive to initial conditions and to display in some 
cases nonmonotonic behavior. Our analysis is centered on cases with /'(0) = 0, and uncovers findings 
qualitatively as well quantitatively different from the more familiar reaction diffusion equations with 
/'(0) > 0. These differences are ascribable to the disparity in time scales between the evolution of 
the front interior and the front tail. 

PACS numbers: 82.40.Ck, 47.20.Ky, 05.45.-a 



I. INTRODUCTION 

Reaction diffusion models are ubiquitous in science. 
They are commonly employed to represent systems 
whose components move diffusively and whose interac- 
tion events, described by the reaction terms, may be rep- 
resented by nonlinear expressions in the macroscopic ob- 
servables such as the system density. Common examples 
can be found in aggregation deposition [2j, chemical 
reactions Q, flame combustion pulse propagation in 
nerves @ and population dynamics [1, 0] ■ Extensions of 
reaction diffusion studies to convective transport [1, Q , 
non-diffusive transport [ToL [Tlj and spatially non-local 
interactions [H, [H| have also been studied in the recent 
past. Here we restrict our attention to a single-species 
reaction diffusion equation in 1-D in its simplest form, 
i.e., 



du(x,t) ^d 2 u(x,t) 



dt 



dx 2 



af(u), 



(1) 



where u(x,t) represents the density profile of the species 
expressed here as a dimensionless quantity, D is the dif- 
fusion constant, a the growth rate and f(u) the nonlin- 
earity. We will further assume that /(0) = /(l) = 0, 
which is a property of the nonlinearity in many systems 
of interest. 

Equations such as (P) often result in propagating wave- 
fronts. The class of reaction terms that allows this fea- 
ture is rather broad but thr ee g eneric types of nonlin- 
earity may be distinguished 14]. One type, henceforth 
called the first type, corresponds to positive f(u) for 
< u < 1 with f(u) w u for u « 0. A well-known ex- 
ample is provided by the Fisher-Kolmogorov-Petrovskii- 
Piskunov (FKPP) equation [H whose reaction term 
is /(it) = u(l — u). Another type, henceforth called the 
second type, corresponds to negative /(it) for < u < b 
and positive f(u) for b < u < 1 such as the Zel'dovich — 
Frank-Kamenetsky [I?} equation (ZF), also referred to 
in the literature as the reduced Nagumo equation [tJ , for 



which f(u) — u{u — 6)(1 — u). This change of sign in 
the nonlinearity is responsible to what is referred to in 
population dynamics as the Allee effect fl8]; a density 
threshold exists below which an initial population even- 
tually gets extinct. Recent work on pattern formation in 
the presence of the Allee effect may be found in ref . [l3| . 
Finally, what we will call the third type of nonlinearity 
has f(u) positive for < u < 1 but is nonlinear in u for 
small u. Reaction diffusion equations with these kinds of 
reaction terms have been used, for example, in studying 
thermal combustion waves [H, [Hi , certain autocatalytic 
chemical reactions (2fjj and calcium deposition in bone 
formation [21j . 

If the reaction term is of the first type, as described 
above, it has been shown that there exists a minimum 
speed for the existence of traveling fronts [TBI, HB] ■ Such 
fronts are termed pulled fronts due to their dynamics be- 
ing dictated by the growth and spreading of the front 
tail [13]. The value of the asymptotic front speed can 
be simply obtained by calculating the spreading of small 
perturbations around the unstable state u = 0. For 
these nonlinearities, initial conditions whose fronts are 
sufficiently steep eventually settle into the traveling front 
shape with speed 2y / Daf'(Q). On the other hand, initial 
conditions with shallower fronts have an infinite set of 
possible front solutions with corresponding speeds [23| . 
These so-called pushed fronts derive their name from the 
fact that the dynamics in the nonlinear region of f(u) 
drives the front propagation [23[ . In such cases, lineariza- 
tion techniques applied to the reaction diffusion equation 
do not allow one to find the traveling front speed. The 
speed selection mechanism differs greatly between the 
pushed and the pulled regimes. The long-time conver- 
gence to the traveling front shape and speed is argued to 
be exponential in the former case (23j and shown to be 
algebraic in the latter (22j . 

Whereas a great deal of attention has been given in 
the reaction diffusion literature to the determination of 
the traveling wave front speed [DljUiJ , the speed selection 
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problem, i.e., the prediction of the asymptotic front speed 
from a given initial condition 25;, 26] , and the long-time 
rate of convergence to that speed [27| , the problem of the 
full transient dynamics has received little attention. Our 
interest is in making some contribution to the study of 
this problem. We ask how an initial condition evolves to 
the asymptotic traveling front shape and speed. We do so 
by analyzing the transient dynamics of reaction diffusion 
equations whose f(u)'s belong to the third type. 

The paper is organized as follows: in Sec. [Til we se- 
lect three specific forms of the nonlinearity and exhibit 
their resultant analytically obtained traveling wave front 
speed and shape. Sec. IIIII is devoted to the study of 
their transient dynamics. Sec. IIVI contains concluding 
remarks. 



II. THREE SPECIFIC NONLINEARITIES AND 
ANALYTIC TRAVELING WAVE SOLUTIONS 



It is known [23, |28| that only pushed fronts exist for 
reaction terms with /'(0) = 0. For some of these nonlin- 
earities one can find exact analytic solutions for the trav- 
eling wave, making possible an exact description of both 
its shape and speed. We exploit the existence of these 
pushed front solutions in studying the transient dynam- 
ics. We consider three examples of positive nonlinearity 



for < u < 1 with f(u) 



0. These are depicted in Fig. [TJ 



as well as f(u) 



for 



0.15 




FIG. 1: Plot of the nonlinearities f(u) under consideration. 
The solid and dashed curves represent f(u) quadratic in u for 
small u as given, respectively, by Eqs. ([2} and ||4j. The dotted 
curve represents f(u) given by Eq. (|6]) with e — 16/(81\/3), 
an example of a nonlinearity cubic in u for small u. 



A. A logarithmic nonlinearity 

Consider what we will henceforth call the logarithmic 
nonlinearity, 



/(«) = ( u + 1) [2 - m(2) + In (u + 1)] 

l)[ln(2)-ln( U +l)], (2) 



In 2 (u - 



which is quadratic in u as u — > and decreases lin- 
early around u = 1. Despite its complicated appear- 
ance, we have managed to find exact analytic solutions 
of the traveling wave front and speed. The procedure [29| 
requires the parametrization of the traveling front equa- 
tion with du/dz instead of z. With the resulting dif- 
ferential equation solved exactly, an analytic value of 
the front speed can then be obtained. Integration of 
(du/dz)^ 1 and its inversion allows the determination of 
the exact tr aveli ng front profile. We have found the speed 
to be c = V Da ln(2) and its traveling front shape (cen- 
tered around the origin) to have the analytic form, with 



u(z) = -1 + 2 



1 i ln(4/3) „z 
T ln(3/2) ^ 



a/D 



(3) 



B. A special case of the ZF nonlinearity 

Another example of a reaction term quadratic in u 
as u — > may be obtained from the ZePdovich-Frank- 
Kamenetsky equation by putting b = 0, i.e., 



f{u) = u 2 {l~u). 
The exact traveling wave front is given by 

u(z) 



1 



z^/a/2D 



(4) 



(5) 



with c = \J Da/2 as the front speed. In the following we 
call this f(u) the quadratic nonlinearity. 



C. A cubic nonlinearity 

A nonlinearity which is cubic in u as u — > can be 
obtained from a variety of systems including, e.g., 



f(u) = esin(7ru) [1 — cos(7ru)] 



(6) 



Given the periodicity of Eq. ([6]), we use it in our study 
only for initial conditions such that u < 1. The speed 
of the traveling front is given by c = v eirDa and its 
traveling front shape is given 



u{z) 



arctan e 



f e ~ z V ena/D\ 



(7) 



The reaction term in Eq. ([6]) is a special case introduced 
to study the dynamics of the angle between the electric 
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field and the polarization in ferroelectric chiral smectic 
liquid crystals [30]. For comparison purposes we choose 
e = 16/(81v3)- This value of e ensures that Eqs. (JT]) and 
(|5|) coincide at their peak u = 2/3. In the following we 
will call the cubic nonlinearity the f{u) in Eq. ([6]) with 
e as given. 



A. Varying the characteristic length 



III. TRANSIENTS DURING THE FORMATION 
OF THE TRAVELING WAVE FRONTS 



Our procedure for the study of transients consists of 

(i) assuming an initial condition roughly in the form of 
a 'right step' with u — 1 as x — > — oo and u = as 
x — > +oo, the form of the switchover being similar in 
shape, but not identical, to the eventual traveling front, 

(ii) computing a quantity that we call the excess speed, 
and (iii) studying the effects of various features of the 
initial conditions on the transient dynamics for the three 
types of nonlinearity considered. The excess speed is the 
amount by which the (eventual) traveling front speed c is 
exceeded by v(t), the instantaneous time rate of change 

of the area under the front, i.e., d J^^dxu(x,t) /dt,as 

explained, e.g., in ref. [8|. It is a convenient quantity to 
focus on because it allows one to calculate how, overall, 
a front profile advances independently of how, on a local 
scale, the profile moves forward. 

Our numerical integration of the reaction diffusion 
equation is performed via an Adams-Bashforth-Moulton 
predictor-corrector method with a spatial mesh of step 
size 0.01 in units of \J D / a, which is the characteris- 
tic length of Eq. ((T|). We consider the front profile as 
separated into a shoulder for most of which u = 1, a 
tail which describes the final vanishing of u at long dis- 
tances, and an in-between interior part. We construct 
our initial conditions by modifying the asymptotic trav- 
eling shape in these three parts separately or in com- 
bination. We report results mainly for tail-modified and 
interior-modified initial conditions since they are the ones 
that show the largest effects. 



In Sec. IIII Al below, we first study initial conditions 
for the front obtained si mply by replacing the character- 
istic length parameter yjD/a in the exact traveling front 
shape by a different value, the front shape being thus 
steeper or shallower than, but qualitatively identical to, 
the exact traveling front. In Sec. IIII Bl and lHI C[ by con- 
trast, we consider initial shapes that are different, even 
qualitatively, from the eventual traveling front shape. We 
obtain them by modifying, in the latter, only the interior 
portion in Sec. IIII Bl and both the interior portion and 
the tail in Sec. IIII CI Finally, in Sec. IIII D| we discuss 
what types of initial fronts will eventually settle into the 
exact traveling front profiles. 



The characteristic length y/D/a of the reaction diffu- 
sion system ((T|) appears naturally in the exact expres- 
sions ©, © and ([7]) for the traveling wave profile for 
the three nonlinearities, and controls the steepness of 
the front. In this subsection we investigate the tran- 
sient occur ring for initial conditions obtained by replac- 
ing \Ja/D by ^^JajD in the respective Eqs. ©, (O 
and ([7]). The quantity we vary is £, the dimensionless 
ratio of the length characteristic of the initial condition 
to yjDja. For values of £ larger (smaller) than 1 , the ini- 
tial spatial profile is shallower (steeper) than the exact 
traveling front. 
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FIG. 2: Long-time decay constant r of the initial front speed 
plotted in units of 1/a for the three nonlinearities depicted in 
Fig. HI The horizontal axis is the dimensionless characteristic 
length £ (see text). The triangles, the circles and the squares 
correspond, respectively, to the logarithmic, quadratic and 
cubic nonlinearity. If £ > 1 (< 1), the initial front profile is 
shallower (steeper) than the asymptotic profile, and the front 
speed decreases (increases) monotonically to the asymptotic 
speed. 

Our finding is that an initial condition constructed in 
this way gives rise to a monotonic decay of the front 
speed to the asymptotic value along with a change of the 
front steepness to the asymptotic shape. The long-time 
behavior appears exponential and can be fitted on a log- 
arithmic scale by a straight line. We calculate in this 
way the characteristic exponent of the decay for each of 
the three nonlinearities. We compare them in Fig. H]by 
plotting the decay time r (reciprocal of the exponent) in 
units of 1/a as a function of the relative steepness ratio 
£. It is evident from Fig. [5] that the decay time for the 
cubic nonlinearity is longer than that for the logarithmic 
nonlinearity. In turn, the logarithmic nonlinearity has a 
decay time longer than that for the quadratic nonlinear- 
ity. This hierarchy in time scales is due to the relative 
strength of the /(u)'s: the larger the value of the non- 
linearity, the faster the rate of change of u. For a large 
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portion of the interval < u < 1, the cubic nonlinear- 
ity has u values smaller than the other two. Similarly, 
the logarithmic nonlinearity has smaller values than the 
quadratic nonlinearity. Moreover, the asymmetry with 
respect to £ = 1 shows that the larger the initial devi- 
ation from the asymptotic shape, the larger the decay 
time [3l|. The fact that shallow initial conditions can 
have larger differences with the asymptotic profile than 
steep ones explains the the increase in the transients in 
Fig. [2] as increases above the value 1. 



B. Modification of the interior of the front: 
nonmonotonic behavior 

Transients dynamics become more complex if the ini- 
tial profile and the traveling front profile are qualitatively 
different. One example studied in this subsection is an 
initial profile obtained by modifying part of the interior 
of the exact traveling shape with a straight line segment. 
This segment starts at the coordinates x = 0, u = 1/2, 
and ends intersecting the exact traveling front profile at 
some point x < (see the insect of Fig. [3j> - In other 
words, the initial profile is shallower than the asymptotic 
one in a limited region of space. Given that the portion of 
the aymptotic profile in the front interior modified to get 
the initial shape is rather small, one might expect a sim- 
ple effect including a monotonic decay as in the previous 
section. Surprisingly, this does not happen. Rather than 
the distortion in the upper part of the interior disappear- 
ing, the shape in the lower part of the interior chnages 
from the exact profile and becomes shallow first. During 
the transient dynamics, this situation corresponds to the 
appearance of a maximum in v(t) as shown in Fig. [3] 
At that instant, the evolution changes and v(t) decays 
monotonically to the asymptotic value as shown in Fig. 
[3J Notice that similar non-monotonic transients with the 
appearance of a maximum in v(t) can be observed if the 
initial profile is taken as a modification of the exact pro- 
file for x > 0, e.g., with an exponential function that 
is shallower than the asymptotic front shape. In other 
words, if the lower, rather than the upper, part of the 
initial front interior is made shallower than the asymp- 
totic shape, qualitatively similar dynamics are observed. 

We define the ratio of the slope of the straight segment 
to the slope of the asymptotic traveling front at u = 1/2 
as the relative steepness a. In Fig. [5J we plot the time 
Ti (in units 1/a) at which the maximum of v{t) appears 
as function of a. The inset shows the value H of the 
relative excess speed maximum as function of a. As a 
approaches 1, the distortion eventually becomes a small 
perturbation. The amplitude of this perturbation decays 
monotonically to the asymptotic profile, similarly to the 
long-time behavior described in Sec. IIII Al This explains 
why, for a > 0.85, the maximum disappears. Indeed, the 
inset shows that the amplitude of the maximum becomes 
smaller as a increases. The other extreme corresponds to 
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FIG. 3: Nonmonotonic behavior of the relative excess speed 
[v(t) — c] /c as function of time. The solid, dashed and dot- 
ted curve represent the evolution of the relative excess speed 
for the logarithmic, quadratic and cubic nonlinearities, re- 
spectively. The initial condition is constructed as follows; it 
is a modification of the interior of the exact traveling front 
profile with a straight line segment depicted in the inset as 
dash-dotted. Only the initial condition for the logarithmic 
nonlinearity is shown in the inset since the three initial con- 
ditions are very similar to each other. 



a approaching when the initial profile tends to flatten 
out for x < 0. In such a case, the dynamics are domi- 
nated by the reaction term since the diffusion becomes 
negligible. It is a simple exercise to show that a flat pro- 
file for u > 1/2 approaches u = 1 exponentially fast. This 
explains why, for a < 0.5, the maximum in the relative 
excess speed also disappears. 
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FIG. 4: Plot of Ti the transient time for the appearance of 
the maximum in the relative excess speed as function of the 
relative steepness a of the initial conditions displayed in the 
inset of Fig. [3] The inset shows the height H of the maximum 
of the relative excess speed as function of a. 
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C. Modification of both the tail and the interior 

A common characteristic of the three nonlinearities we 
study is /'(0) = 0. This implies that the dynamics in 
the tail of the front are slower than the dynamics in the 
front interior. Therefore, we may expect an even more 
complex transient if the tail as well as the interior of the 
front are modified. We show our numerical results for 
this situation in Fig. [5] for the logarithmic nonlinearity. 
The initial condition is chosen as follows. The lower part 
of the interior of the exact traveling front is cut off at 
an amplitude A = 0.1 by a straight line whose length 
depends on its steepness. This cut-off thus creates a tail 
which is infinitely steep. At the same time the rest of 




t (units of 1/a) 



FIG. 5: Relative excess speed in the case of the nonlinearity 
given by Eq. (J2| for an initial condition constructed as ex- 
plained in the text. The various insets shows how the evolving 
front is different from the exact traveling front shape (dotted 
line) at different times. The exact traveling front is plotted 
by making the two curves cross each other at u — 1/2. 

the front interior is modified from the shape of the exact 
solution by increasing the characteristic length yj D/a as 
in Sec. IIII Al In other words, the initial condition has 
an interior part which is shallower than the exact trav- 
eling front with its lower part modified with a straight 
line. Such an initial condition allows one to observe three 
different transient regimes. The first is associated with 
the time necessary for the interior part of the front to re- 
duce its steepness and become essentially indistinguish- 
able from the traveling front shape (see the center inset 
in Fig. [5] The second transient is created by the fast 
upwards propagation of the shallow straight line profile 
(see the top left inset in Fig. [5]) present initially in the 
lower part of the front interior. These dynamics are sim- 
ilar to the one observed in Figs. [3]and|4j eventually the 
interior front shape is converted into a shallow profile. 
When that happens once again v(t) reaches a maximum. 
The final transient regime then corresponds to the mono- 
tonic decay to the asymptotic speed as shown earlier in 
Fig. [2 Notice that, even though the initial front has a 
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t (units of 1/a) 

FIG. 6: Plot of the relative excess speed for the logarithmic 
nonlinearity as function of time for four initial fronts with 
different value of A (see text). The solid, dotted, dashed and 
dash-dotted curves correspond to A — 0.04, 0.06, 0.1 and 0.12, 
respectively. The inset shows how the delay time T2 between 
the minimum and the maximum changes as function of A for 
the three nonlinearities. 



tail with compact support, i.e., is infinitely steep, that 
extreme steepness is not transferred to the front interior. 
On the contrary, it is the shallow interior that eventually 
makes the interior profile shallower than the traveling 
wave shape during the transient. 

The appearance of the two extrema is also studied by 
selecting initial conditions whose interior parts have a 
given characteristic length larger than the asymptotic 
traveling front, and whose right-most portion is made 
up of a straight line segment proceeding from u = A to 
u = 0. We keep the projection of the straight segment on 
the x-axis fixed and vary A. In the main part of Fig. [51 
we display four curves (solid, dotted, dashed and dash- 
dotted) with, respectively, increasingly larger values of 
A. The effect of the time scale disparity between the tail 
and front dynamics can be appreciated from Fig. [5] At 
first, as A increases, the time lag between the minimum 
and maximum also increases. However, beyond a certain 
value of A, a larger portion of the shallow profile sees a 
rapidly increasing shape of the nonlinearity f(u). The 
evolution of the initial front thus gets increasingly faster 
as A increases, reducing the time necessary for the profile 
to become shallow overall. This explains why the maxi- 
mum of the dashed and dash-dotted curves move to the 
left. The inset of Fig. [5] shows the variation of this time 
lag called T2 as function of A for the three nonlinearities. 
This transition is visible in all three cases but is more 
pronounced in the logarithmic and quadratic nonlineari- 
ties. 
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D. Initial fronts with long tails 

Further studies of the transient dynamics we have car- 
ried out have shown that the types of initial fronts that 
eventually settle into the traveling wave shape is quite 
large. For the three nonlinearities considered, initial 
fronts whose tails decay as a power law with exponent 
> 3, approach the asymptotic speed c within the accu- 
racy of our numerical scheme [32|. For shallower initial 
tails, the asymptotic speed seems to settle to larger val- 
ues indicating that the basin of attraction of the traveling 
front solution not only contains initial fronts whose tail is 
steeper but also those that are much shallower than the 
asymptotic ones. This is different from the class of non- 
linearities with /'(0) > in the pushed regime. In those 
cases, if the initial tail decays faster than an exponential 
whose characteristic length depends on f(u), the front 
speed converges to the lowest (pushed) asymptotic trav- 
eling speed. On the other hand, any initial conditions 
whose tail is shallower than the above exponential will 
converge to a larger speed value [29j . This is not the case 
in the nonlinearities we have studied since the front inte- 
rior dynamics is faster than the tail dynamics. The initial 
shape of the front interior is thus more important than 
the initial tail shape. Our various examples have shown 
that the initial front adjusts itself to the traveling front 
depending mostly on the shape of the interior part of the 
front. The initial tail seems to determine the asymptotic 
speed only if it vanishes at long distances with a shallow 
power law profile. 

IV. CONCLUDING REMARKS 

We have studied the transient dynamics of reaction dif- 
fusion systems whose reaction term f(u) is characterized 
by being positive for < u < 1 with /'(0) = 0. Specif- 
ically, we have chosen two nonlinearities whose f{u) is 
quadratic and one whose f(u) is cubic for u ps 0. We have 
found analytic expressions for the corresponding travel- 
ing front shape and speed. With the help of these expres- 
sions, we have analyzed how certain types of initial con- 
ditions evolve in time. Our initial conditions are simple 
modifications of the exact traveling front shape. Mono- 
tonic relaxation to the traveling front speed is observed if 
the initial front is qualitatively similar to the asymptotic 
profile, being constructed from the latter only by varying 
the characteristic length scale of the traveling wave front. 
Nonmonotonic behavior instead may be observed when 
the initial conditions are qualitatively different from the 
asymptotic profile. In such cases, the transient dynamics 
is dominated by the evolution of the interior part of the 
front since the tail evolution is much slower due to the 



nonlinear nature of f(u) for u — > 0. If part of the inte- 
rior is shallower (steeper) than the asymptotic traveling 
front, the shallowness (steepness) is transferred to the 
whole interior profile, initially. Subsequently the speed 
then decays monotonically to zero with the long time 
behavior being exponential. 

Many reaction diffusion systems, including the well- 
known FKPP case [ID, HH, have /'(0) ^ 0, which means 
that, if u is the density of a species, growth can occur 
without the interaction of members of the species. On 
the other hand, our analysis in this paper has centered 
on cases in which /'(0) = 0, which represents situations 
in which even at low concentration, interaction (nonlin- 
ear term in u) is essential to growth. We find major 
differences relative to systems with /'(0) > 0. In these 
latter systems, in the pushed regime, the initial fronts 
that eventually settle into the asymptotic shape are gen- 
erally those whose tail decays in space exponentially with 
a characteristic length smaller than a value dependent on 
the form of f(u). In our investigation, on the other hand, 
the initial conditions that settle into the traveling front 
can have tails as shallow as power laws. Given the time 
scale disparity between the dynamics of the tail and of 
the interior, the tail evolution can be said to be slaved 
to the interior evolution. For this reason, the tail does 
not affect the speed selection mechanism unless it is ex- 
tremely shallow. As the integer n increases in reaction 
terms with f(u) ss u n close to u = 0, the disparity in 
time scales grows bigger and we expect the basin of at- 
traction for the initial fronts to get larger. With the 
reaction terms that we have used it is difficult to verify 
numerically that this is indeed the case. However, we 
can easily make a comparison with a nonlinearity whose 
tail dynamics are even slower than in three cases we have 
studied. By taking the ZF nonlinearity with b ^ 0, we 
clearly see that initial conditions that decay with power 
laws of exponent > 2 approach the traveling front profile. 

In this work we have only scratched the surface of the 
important topic of transient dynamics. We hope that our 
study will motivate further research on the formation of 
the traveling front profile in reaction diffusion systems. 
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